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ABSTRACT 

In a recent paper we have studied the redshift power spectrum P s (k,fi) in 
three CDM models with the help of high resolution simulations. Here we apply 
the method to the largest available redshift survey, the Las Campanas Redshift 
Survey (LCRS). The basic model is to express P s (k,fi) as a product of three 
factors 

P s (k^) = P R (k)(l + p^) 2 D(k^). (1) 

Here \i is the cosine of the angle between the wave vector and the line of 
sight. The damping function D for the range of scales accessible to an accurate 
analysis of the LCRS is well approximated by the Lorentz factor 

£) = [l + l(fc/i(7 12 ) 2 ]- 1 - (2) 

We have investigated different values for (3 ((3 = 0.4, 0.5, 0.6), and measured the 
real space power spectrum P R {k) and the pairwise velocity dispersion au(k) 
from P s (k,fi) for different values of \i. The velocity dispersion a 12(h) is nearly 
a constant from h = 0.5 to 3 /iMpc -1 . The average value for this range is 
510 ± 70kms _1 . The power spectrum P R (h) decreases with h approximately 
with k~ 1,7 for h between 0.1 and 4 /iMpc^ 1 . The statistical significance of the 
results, and the error bars, are found with the help of mock samples constructed 
from a large set of high resolution simulations. A flat, low-density (Qq = 0.2) 
CDM model can give a good fit to the data, if a scale-dependent special bias 
scheme is used which we have called the cluster-under-weighted bias (Jing et 
al.). 



Subject headings: galaxies: clustering - galaxies: distances and redshifts - 
large-scale structure of Universe - cosmology: theory - dark matter 
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1. Introduction 



The spatial distribution of galaxies can be retrieved approximately, when the redshift 
is used as a measure of the distance. It is not exact, because the galaxies in general deviate 
from the linear Hubble flow, but the redshift distortions can provide us with valuable 
information about the dynamics of galaxies. A traditional method for studying the redshift 
distortion is the redshift two-point correlation function (RTPCF) £ 2 (r p , tt) flGeller fc Peebles 
15751 , Peebles EM], |Bean et al. 1553] Pavis k Peebles 19831 |Mo, Jing fc Bonier 1993 
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fc Gcllcr 1999| . [Small et al. 1999 ). Assuming certain functional forms for the distribution 
function (DF) of the pairwise velocity (say, an exponential form, Peebles 1976| ) and for the 
average infall velocity (say, the form of the self-similar solution, Pavis fc Peebles 1983|) , one 
can construct a model for t; z (r p ,ir) which describes the real situation, when the coupling 
between the peculiar velocity and the spatial density of the galaxies is weak (Peebles 1980| ) . 
A comparison of the model with observations for £ z (r p ,7r) provides a test for the validity 
of the model assumptions (such as the DF of the pairwise velocity) and a determination 
of the pairwise velocity dispersion aia(r). Early studies of small redshift samples ( Pccblcfj 
|1976[[P"eebles 1979| ) have shown that an exponential form for the DF of the pairwise velocity 
is preferable to a Gaussian form. This important conclusion has been confirmed by a 
number of later studies based on much larger surveys ( Pavis fc Peebles 1983[ [Marzke et 
al. 1995| , |Fisher, et al. 1994] ). The pairwise velocity dispersion was not easy to pin down; 
the value sensitively depends on the presence or absence of rare rich clusters in the redshift 
surveys which contained only a few thousand galaxies ( [Mo, Jing fc Borner 1993] , [Zurck ct| 
|al 1994[ [Marzke et al. 19951 ). Motivated by this fact, several authors have since attempted 
( [Kepner, Summers fc Strauss 19971 pavis, Miller fc White 19971 , IStrauss, Ostriker fc Cenl 
|1998[ [Baker, Davis fc Lin 2000] ) to find new statistics for the thermal motion of galaxies 
which are less sensitive to the regions of rich clusters. The new statistics, by their design, 
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are more robust with respect to the sampling of rich clusters than the pair- weighted 012 (r) 
and have produced interesting constraints on models of galaxy formation, but the results 
sometimes are more difficult to interpret in the context of dynamical theories (e.g. the 
Cosmic Virial Theorem and the Cosmic Energy Equation). Fortunately, the pair- weighting 
will no longer be a problem for future large redshift surveys like 2dF and SDSS which will 
contain hundreds to thousands of rich clusters in the survey. In fact, the currently largest 
publicly available redshift survey, the Las Campanas Redshift Survey (LCRS; |Shectma"n 



|et al. 1996 ) which contains about 30 rich clusters, is already large enough to reduce the 



sampling error in <J\2{t) to 15% only ( |Jing, Mo fc Borner 1998| , hereafter JMB98). Therefore 
the pairwise velocity dispersion, which is a well-defined physical quantity and has simple 
relations to the dynamical theories (the Virial Theorem), will remain an important quantity 
to measure in observations. 

Although the pair- weighting is not a problem for the statistics of 0"i2(?") with the largest 
redshift surveys or the next generation of redshift surveys, the currently widely used method 
of measuring a^ir) QDavis fc Peebles 1983|) does have its limitations (JMB98). Functional 
forms must be first assumed for the DF of the pairwise velocity and the infall velocity, and it 
must be assumed that the spatial density and peculiar velocity of galaxies are uncorrelated 
so that the redshift two-point correlation function (model) can be written as a convolution 
of the real-space correlation with the DF of the pairwise velocity. Although the validity 
of the functional form for the DF of the pairwise velocity can be checked by comparing 
the model of the RTPCF to the observed one, the tests of different functional forms are 
rather limited. The functional form may depend on the separation between the objects 
( |Juszkiewicz, Fisher fc Szapudi 1998| , [Magira, Jing fc Suto 2000]) . The infall velocity also 



has a significant effect on the determination of cri2(r), because the two quantities are rather 
degenerate in modeling £, z (r p , n) flJing fc Borner 1998|) . Moreover, the density and velocity 
are strongly coupled at the small scales where previous studies have measured a^ir). These 
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limitations can lead to some systematic bias in the estimation of 012(7") (JMB98; Jing fc 



Borner 1998| ). These systematic effects were not so important in early statistical studies 
of 012(7") since the sampling errors in 012(7") were dominant. Since the sampling error is 
just 15% for the LCRS (JMB98) and will be even smaller for the upcoming surveys, it is 
important to consider alternative methods to measure the redshift distortion and 012(7") 
which are less model dependent. 

In a recent paper (|Jing fc Borner 200 1| , hereafter JB2001), we have studied the 
redshift power spectrum P s (k,fi) in three typical CDM cosmological models, where \x is 
the cosine of the angle between the wave vector and the line-of- sight. Two distinctive 
biased tracers, the primordial density peaks ( Bardcen et al. 1986 ; Davis et al. 1985 ) and 



the cluster-under- weighted population of particles (JMB98), are considered in addition 
to the pure dark matter models. Based on a large set of high-resolution simulations, we 
have measured the redshift power spectrum for these three tracers from the linear to the 
non-linear regimes. For all the tracers and the three CDM models, the redshift power 
spectrum P s (k,fi) can be expressed by 

P S (k,fi) =P R (k)(l + (3fi 2 ) 2 D(k f icr 12 (k)) (3) 

where P R (k) is the real space power spectrum, and (3 = Q os /b (Q is the density parameter 
and b is the linear bias parameter). The factor (1 +/3/z 2 ) 2 accounts for the linear compression 
of structures on large scales ( |Kaiser 1987| ). Although 012 (&) and the 3-D peculiar velocity 
dispersion are not immediately comparable, it was pointed out in JB2001 that their values 
are different only by 15% in the simulation, if r = 1/k is used in the comparison. So 012 (k) 
is a good indicator for the velocity dispersion. In Equation ([3D the damping function D, 
which should generally depend on k, \x and 012 (A;), is a function of one variable k^ia^k) 
only. The functional form of D(k/j,ai2(k)) was found to depend on the cosmological model 
and the bias recipes, but for small k (large scales) where D > 0.1, Dik^ia^k)) can be 
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accurately expressed by the Lorentz form, 

D(k,a 12 (k)) = i - W ^ W) ■ (4) 

This implies that the exponential form is a good model for the distribution function of the 
pairwise velocity at large separation. Applying the model of Equation (|3|) to an observation 
of P s (k,fi) therefore may yield a determination for the four quantities: P R (k), (3, a±2(k) 
and the functional form of D(y) all of which are useful observables for testing galaxy 
formation models. An additional advantage is that the infall effect is accounted for by the 
single parameter (3, compared with a complicated form for the infall velocity used in the 
correlation function method. 

In this paper, we apply this method to the LCRS, the largest redshift survey publicly 
available. The redshift power spectrum P s (k,fi) is measured from the RTPCF £ z (r p ,7r) of 
this survey (JMB98) using the method outlined in Section 2. The result is presented in 
section 3. The statistical significance of the measurement, and the error bars, is tested using 
mock samples in section 4. The final section (§5) is devoted to discussion and conclusions. 

2. Method 

We convert the redshift two-point correlation function £ z (s) to the redshift power 
spectrum by the Fourier transformation: 

P 5 (k) = |^(s)e- tk - s cfe. (5) 

In cylindrical polar coordinates (r p , <ft, it) with the 7r-axis parallel to the line-of-sight, 
P 5 (k) depends on k p , the wavenumber perpendicular to the line-of-sight, and on k n , the 
wavenumber parallel to the line-of-sight. The power spectrum can be written 

P s {k p ,k n )= h z {r p ^)e- i[k ^ cos ^ )+ ^\- v dr p d<pdTi . (6) 
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With some elementary mathematical manipulation, we get the following expression: 

/oo /*oo 
dir \ r p dr p £ z (r p , tt) cos(k w ir) J (k p r p ) (7) 
-oo JO 

where Jo{k p r p ) is the zeroth-order Bessel function. 

We will use the redshift two-point correlation function £ 2 (r p ,7r) of the LCRS measured 
by JMB98. In their paper, JMB98 analyzed the correlation function for the whole sample as 
well as for the three southern slices and the three northern slices separately. Moreover, they 
have generated a large sample of mock surveys with different observational selection effects 
and measured £ z (r p , tt) for them. Their tests with the mock samples are very important, as 
emphasized by JMB98 themselves, for checking the statistical methods, correcting for the 
observational selection effects, and estimating the statistical errors of the measurements. 
We will also use their results of the mock samples in this paper for the same purposes. 

In their work, JMB98 measured £ 2 (r p ,7r) in equal logarithmic bins of r p and in equal 
linear bins of tt. The reason why different types of bins are chosen for r p and tt is the 
fact that £z(V p ,7r) decreases rapidly with r p but is flat with tt on small scales. Thus this 
way of presenting ^ z {t p ,tt) is better than using the log-log or the linear-linear bins for r p 
and tt, and is also suitable for the present work. The peculiar velocity of a few hundred 
kms -1 should smoothen out structures on a few /i _1 Mpc in the radial direction, and the 
linear bin of Att^ = 1 h~ 1 Mpc is suitable for resolving the structures in the radial direction. 
With logarithmic bins chosen for r p , the r p dependence is resolved well, because otherwise 
the small scale clustering on the projected direction cannot be recovered. With this bin 
method, we obtain the power spectrum: 

P s (k p , k n ) = 2tt A7r i r pj A mr P,j£z( r P,j> ft) cos^TT^Joikprpj) (8) 

where 7Tj runs from —50 to 50/i _1 Mpc with Atti = 1 /i -1 Mpc and r p j from 0.1 to 
31.6 /i _1 Mpc with Alnr p j = 0.288 (Be careful not to confuse two its in Eqs. (0) and (|8|): 
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the first 7r in the right-hand-side has the conventional meaning, i.e. 3.14159..., and the 
others are for the axis along the line-of-sight.). We use the summation of Eq.(8) with 
rectangular boundaries in n and r p . This gives good, unbiased results for both the velocity 
dispersion and the power spectrum, as will be seen in §4. This works also better than a 
spherical window of radially decreasing weighting which tends to bias the power spectrum 
estimates. We have tested this throughly comparing the results for mock samples with 
those for the full simulation. 

There are several sources which could introduce errors to our measurement of the 
redshift power spectrum. One is the Poisson error related to the discrete sampling of 
galaxies, and the other is the cosmic variance related to the limited volume of the survey. 
Both types of uncertainty are intrinsic for any redshift survey, and are already factored into 
the measured £ z (r p ,7r). When converting the correlation function to the power spectrum 
with Eq.(H), the finite bins and the cutoffs of r p and 7r may cause additional errors. To 
reduce the effect of the finite bins, we divide every 7r and lnr p bin into N sub-bins, and 
use the cubic spline method to interpolate £ z (r p ,7r). We have run a number of trials for 
different numbers of the sub-bins N, and found that the final result is insensitive to N. For 
the results presented in the following sections, N = 10 is used. Since these different types 
of errors are mixed in a complicated way in the measured power spectrum, we will use the 
mock samples to quantify the errors. 

Landy, Szalay & Broadhurst (1998) have also considered the Fourier transform of the 
redshift two-point correlation function for the LCRS survey. Although the spirit of their 
work and ours is similar, there are some important differences. We have included the infall 
effect in our statistics and consider the full angular /i-dependence of the redshift power 
spectrum. Furthermore, we will use the two dimensional integral (Eq.[7]) to determine both 
the velocity dispersion and the power spectrum. 
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3. The redshift power spectrum in the Las Campanas redshift survey 

Since the redshift distortion does not change the total number of the galaxy pairs at 
a certain projected distance, it is not difficult to prove from Eq.(3) that P s (k, 0) = P R (k). 
Thus the damping factor D(k,n) can be estimated by the following expression 

The quantities on the right hand side can be measured through Eq.(|8|) except for the 
parameter (3. If the observational catalog is large enough, one can simultaneously determine 
the damping factor and (3 as described by JB2001. However, the LCRS is too small 
to measure this quantity because, as we can see in Figure 1, there is a large statistical 
fluctuation on large scales (ka v fi < 1). This is consistent with the recent estimates of 
Matsubara et al. ( pOOOf) based on an eigenmode method, who give (3 = 0.30 ± 0.39. 
Therefore, we will not attempt to determine the value of (3. Instead we will focus on 
measuring the damping function. Here we fix a value for (3, (3 = 0.5, which is consistent 
with the observations of the cluster abundance (e.g. [Kitayama fc Suto 1996|) as well as with 



the statistical result of Matsubara et al (|2000| ) for the LCRS. We have checked the results 
for (3 = 0.4 and (3 = 0.6, and found that our conclusions are little changed if (3 varies within 
this range. 

Figure 1 shows the inverse of the damping function D{k,fi) we measured for the Las 
Campanas redshift survey. Following JB2001, the damping factor is expressed as a function 
of the scaling variable k[ia v . For this figure, a v is fixed to be a constant 500kms _1 , and 
only the data points which have a relative error smaller than 50% are presented. The 
error of D(k,fi) is estimated from the scatter of this quantity between the cluster-weighted 
mock samples (see §4 about the mock samples). From the figure, it is easy to see that our 
measured result of D(k,/j,) can be well described by a Lorentz form (the solid line). This 
result is consistent with the model studies of JB2001 who have shown that the damping 
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factors of different bias models are all very close to the Lorentz form when the damping 
factor D(k,fi) > 0.1. The LCRS is still not large enough to explore the highly non-linear 
regime where D(k,fi) <C 0.1. This also means that an exponential distribution function for 
the pairwise velocities is a good approximation in this range. In addition, we note that 
there exists significant scatter on large scales k/j,o~ v < 1, thus it is very difficult to measure 
the f3 value accurately with this catalog (cf. Matsubara et al. |2000| ). 



As Figure 1 shows, the damping function that the LCRS can explore is in the interval 
0.1 < D(k,fi) < 1. The accuracy is limited by the sample-to-sample error (the cosmic 
variance) on large scale (D(k,fi) ~ 1). Since the LCRS measured the redshift for each 
1.5 x 1.5 deg 2 field by one exposure, the dense regions in the galaxy distribution (e.g. 
clusters) are under sampled, which seriously limits the accuracy of our measurement for 
small scales [D{k,fi) <C 1]. For the regime which the LCRS can effectively explore, the 
Lorentz form is a good approximation for the damping function, consistent with the model 
study of JB2001. Thus, we can determine the real-space power spectrum P R {k) and the 
pairwise velocity dispersion a v by a least-square fitting to the observed redshift power 
spectrum, assuming the Lorentz form for the damping factor. Although P R (k) = P s (k, 0) 
mathematically, for a finite-sized survey like the LCRS, P (k, 0) fluctuates around the true 
value of P R (k) for the limited number of the modes. We can better determine P R (k) if we 
combine P s (k,fi) at different angles /i, thus we treat P R (k) as a free parameter. 

We show such best fitting curves (the solid lines) together with the measured redshift 
power spectrum in Figure 2. For the curves from top to bottom, the wavenumber k is 
incremented by A\gk = 0.2 from 0.25 /iMpc -1 to 2.5/iMpc -1 . As the figure shows, the 
model we use for P s (k,[i) (i.e. the Lorentz form for the damping function and equation |3|) 
succeeds in describing the result of the LCRS data. Our best fitting values for au(k) and 
P R (k) are plotted in Figures 3 and 4 respectively. Here we have corrected the fiber collision 
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effect of the survey QShectman et al. 1996|) for both quantities following the procedure 



of JMB98. The error bars shown in these figures are simply the standard deviations for 
identical analyses of the 60 mock samples. We have also used the bootstrap method to 
estimate the errors, and found that the bootstrap errors are typically 30 ~ 50% smaller than 
the errors determined from the mock samples. We prefer to adopt the mock errors, because 
they have adequately included the sample-to-sample variations (i.e. cosmic variances). 
In addition, the mock errors are not sensitive to the bias model used; the mock samples 
constructed from the dark matter particles give similar error estimates (cf. JMB98). The 
velocity dispersion 012 (A;) is nearly a constant from 0.5 to 3/iMpc -1 . The average of 012 (A;) 
for this k range is 510 ± 70kms _1 . This error bar for the averaged 012 (A;) is only slightly 
smaller than the typical error ~ gOkms" 1 of each individual k bin, since the errors of 
individual k bins are correlated. Our result for 012 (A;) is in good agreement with our 
previous determination 570 ± 80kms _1 at a separation of r = 1 /i _1 Mpc based on a fit to 
the two-point correlation function (JBM98). A reasonable change for the (3 value from 0.4 
to 0.6 based on the cluster abundance observations (e.g. Kitayama fc Suto 1996| ) leads to a 



change of 5 ~ 10% in 012 (A;). The south (dashed lines) and north (dotted lines) subsamples 
give a result consistent with the whole sample within the statistical uncertainty. 

Our measurement for the power spectrum P R (k) is presented in Figure 4. The 
spectrum is approximately a power law for the range of scales that the LCRS can explore. 
It decreases with k approximately as oc A; -1 ' 7 for 0.1 < k < 4/iMpc -1 . We consider this 
result very reliable, because we have tested our method of a direct Fourier transform 
extensively with mock samples. In our tests we found that this method is superior to using 
other types of window functions, e.g. spherical ones, especially when the aim is to obtain 
unbiased estimates of 012 and P{k). These results are qualitatively in agreement with the 
real space power spectrum of the APM survey flBaugh fc Efstathiou, 1994]) . The shoulder 
at k « 0.4/iMpc -1 apparent in Figure 4 appeared also in the APM results (their Figure 11) 
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at exactly the same scale, though the error bars of the LCRS are much larger. We also note 
that the power-law index of P(k) for 0.4 <k< 4/iMpc -1 is » —1.2 compared with » —1.3 
of the APM survey. 

4. Mock samples and test of the method 

As we discussed in Section 2, there are several possible sources which could introduce 
errors to our measured redshift power spectrum and to our fitted quantities. The use of 
the finite bins and the cutoffs in the radial and projected separations may lead to a certain 
systematic bias in the measured power spectrum. When we write down equation (0), we 
implicitly make the distant observer assumption. There are also statistical uncertainties in 
the power spectrum which are related to the cosmic variance and the Poisson error inherent 
to a redshift survey. All these uncertainties, systematic or statistical, are however difficult 
to model with an analytical method. Fortunately, they can be easily quantified with the 
help of the mock samples. 

We will use the mock samples of JMB98 which were produced from a set of CDM 
N-body simulations of 2 million particles. Three flat CDM models were considered: one 
is the once standard CDM model of the density parameter Oo = 1 and the normalization 
erg = 0.62, and the other two are low-density flat CDM universes with parameters (f2 , erg) = 
(0.2, 1) and (0.3, 1) respectively. The known selection effects of the survey, including the sky 
boundaries, the limiting magnitudes and the redshift selection function, the observed rate 
due to the limited number of fibers, and the missing galaxies due to the fiber mechanical 
collision, were properly simulated in the mock samples. Because the CDM models were 
found to have too steep a two-point correlation and too large a pairwise velocity dispersion 
compared to the LCRS result, JMB98 have introduced a simple bias model which under 
weights the cluster regions in the dark matter models. The mock samples based on this 
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bias model yielded much better matches to the LCRS observation, and the model of 
(f2 ,cr 8 ) = (0.2, 1) was shown to be consistent with the observations. 

To test our statistical method, we will use two types of mock samples of the cluster 
weighted Q = 0.2 sample. To test if there is any systematic bias with the method, we will 
use mock samples which incorporate all selection effects except the fiber collision. The 
reason why we exclude the fiber collision effect is because it is not possible to include this 
effect in a 3-D simulation sample, as the results of the 3-D sample must be compared with 
those of the mock samples in the test. To quantify the statistical errors of our measured 
quantities, we will use the mock samples based on the cluster-weighted Qq = 0.2 sample 
with the fiber collision effect. For these mock samples, we will include or exclude the fiber 
collision effect in order to measure the systematic effect of the fiber collision. 

Figure 5 presents our determination of the damping function in the cluster-weighted 
mock samples of Qq = 0.2 according to Equation (0). The right panel shows the result for 
one mock sample which is randomly chosen. The damping function can be approximated 
by a Lorentz function (the solid line) within the scatter as in the LCRS observation. 
If we increase the number of mock samples to 60 and thus reduce the noise, the mean 
damping function decreases slightly faster than the Lorentz form at kfia v > 3 (left panel). 
This behavior was also found by JB2001 for the same model who measured the damping 
function by FFT for the full simulations (thus achieving higher accuracy). In fact, the 
mean damping function is not only qualitatively but also quantitatively in good agreement 
with the results of JB2001. This test indicates that our determinations for the damping 
function are not systematically biased, though there still exist considerable scatter for the 
LCRS. With the upcoming large surveys 2dF and SDSS, we expect that the uncertainties 
of the redshift power spectrum will be considerably reduced and it will become practical to 
measure a damping function which can be used to discriminate between theoretical models. 
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For the LCRS, however, the Lorentz form is still sufficiently accurate for a description of 
the damping function. 

Using the Lorentz form for the damping function, we determine the pairwise velocity 
dispersion a 12 (k) and the real-space power spectrum P R (k) for the mock samples of the 
Qq = 0.2 model in the same way as for the real catalog LCRS. These two quantities are 
measured for each of the 60 samples, and the mean values and the expected la errors of 
the mean values are plotted in Figure 6 (symbols with error bars). The agreement is quite 
satisfactory with the same two quantities measured from the true particle velocity and 
positions in the simulation (solid lines, k — 1/r for the pairwise velocity dispersion). The 
plots in Figure 6 demonstrate that this method to estimate ai 2 (k) and P R (k) works well, 
and does not suffer from a strong bias. There does exist some small systematic bias, as 
appeared in the estimated velocity dispersion a± 2 (k) at A; ~ 0.3 hMpc' 1 and k « 1.5 /iMpc -1 , 
which is likely caused by the sharp cutoffs in the radial and projected separations and can 
be removed by estimating RTPCF to larger separations. 

Next we compare our statistical results with the predictions of the CDM models using 
the cluster- weighted bias. The results are shown in Figure 7. For the model with fi = 1, 
we have shifted the power spectrum by the factor 1/cxf, assuming a linear bias of l/a$. The 
power spectra are consistent with the LCRS observation for the range of scales explored 
by the survey, but the velocity dispersions of the CDM models are still systematically 
higher than the observation. Only the CDM model with the density parameter Q = 0.2 
is consistent with the observations. Its la lower limit actually coincides with the LCRS 
result. Just as in JMB98 we may conclude that the flat, low-density models give a better 
fit to the data than high-density CDM models. The cluster-weighted bias improves the 
agreement between models and data remarkably. From Fig. 6 we can see that without 
the cluster-weighted bias (dotted lines) the velocity dispersion would increase by roughly 
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100 km s -1 , and the power spectrum in real space would be about 2 times higher. Similar 
changes could be seen in Fig. 7, if the cluster- weighted bias were removed from the model 
calculations. 

In a previous paper (JB2001) we have shown that the bias model using primordial 
density peaks gives even higher values for a^ik) and P R (k), since the clusters have more 
weight than the dark matter particles in that scheme. These arguments carry over to the 
analysis of the LCRS data. It seems that a preference for the cluster-weighted bias emerges 
from our analysis. 

5. Discussions and conclusions 

Here we want to comment briefly on some previous work on the determination of the 
pairwise velocity dispersion <t 12 for the LCRS. In JMB98 a u has been determined based 
on the traditional method of fitting the RTPCF. A value of a 12 = 570 ± 80 km s _1 at 
r p = 1 h~ l Mpc has been obtained, but figure 1 of JMB98 shows that this quantity is slightly 
smaller (about 450 to 500 kms -1 ) at smaller (r p < 0.3/i _1 Mpc) or larger (r p > 3 /i _1 Mpc) 
separations. Since some difference between the quantities in Fourier space and coordinate 
spaces is expected, we regard our result in this paper, <r 12 = 510 ± 70kms~ 1 around 
k = 1/iMpc -1 , to be in excellent agreement with the result of JMB98 as well as the 
recent results of Ratcliffe et al. (1998). The value is larger than that Landy et al (1998) 
obtained based on an power spectrum analysis of the LCRS who ignored the infall effect 
and considered P s (k,n) at fi = (0, 1) only. In addition, the estimate of the power spectrum 
P R {k) for the LCRS as a direct derivation from the redshift power spectrum seems a reliable 
result in the interval between 0.1 /iMpc -1 and 3 /iMpc^ 1 . We suspect that guesses at P R (k) 
for larger scales (smaller k) suffer from the sample variation of the LCRS. 
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We reach the following conclusions: 

• The power spectrum P R (k), and the pairwise velocity dispersion cr^ik) can be 
estimated reliably from the redshift power spectrum of the LCRS. An assumed form 
for the distribution function of the pairwise velocity is not needed, although it turns 
out that the damping function fits reasonably well to the Lorentz form, over an interval 
in k from 0.1 /iMpc -1 to 3/iMpc -1 . The Lorentz form for the damping function is the 
Fourier transform of the exponential distribution of the pair velocities. The fact that 
over a range of the scales the Lorentz form is a satisfactory approximation shows that 
in a corresponding range of r-scales, the exponential distribution is a good model. 

A value for [3 must be assumed, but it has been reasonably determined by other 
observations. 

• The last values are a 12 (k) = 510 ± 70kms _1 from k = 0.5 to 3 /iMpc^ 1 , in agreement 
with JMB98; P R (k) oc fc" 1 - 7 for k = 0.1 to 4/iMpc -1 . 

• Very important is the extensive use of mock samples which are constructed from a 
large set of high resolution simulations. They allow us to estimate the errors and the 
statistical significance of the results reliably, despite the intricate way which various 
error producing effects interact. 

• We find reasonable fits to these measured quantities only in flat, low-density CDM 
models with a cluster-weighted bias. This antibias model-the number of galaxies 
per unit mass in a massive cluster is proportional to M a (M is the cluster mass) - 
suggested in JMB98 obviously deserves attention in future work. It should be pointed 
out that the parameter a = —0.08 used in this paper corresponds to a particular 
selection of galaxies: the LCRS galaxies. For other mix of galaxies (i.e., other 
observational criteria), a could be different. In our recent work (in preparation), we 
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found the bias model of a — —0.25 can successfully account for the clustering of IRAS 
galaxies. 

• The method applied here would be more powerful in samples larger than the LCRS. 
The power spectrum P R {k), the /3-value, pairwise velocity dispersion, and the 
damping function can all be determined simultaneously, if the data are good enough. 

We are grateful to Yasushi Suto for the hospitality extended to us at the physics 
department of Tokyo university where most of the computation was completed. J.Y.P. 
gratefully acknowledges the receipt of a NAO COE foreign research fellowship. The work is 
supported in part by the One-Hundred- Talent Program, by NKBRSF(G19990754) and by 
NSFC to Y.P.J., and by SFB375 to G.B.. 
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Fig. 1. — The inverse of the damping function D(k, //) of the Las Campanas Redshift Survey, 
plotted as a function of k\xo v . In this plot, we have taken (5 = 0.5 and a v = 500 km s -1 . 
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Fig. 2. — The redshift power spectrum P s (k,fi) of the Las Campanas Redshift Survey. The 
solid lines are the best-fit curves for free parameters a v and P R (k), assuming the Lorentz 
form for the damping function. Different colors correspond to different values of k. From 
top to bottom, k = 0.25 to 2.5/iMpc -1 at A lg ft = 0.2. 
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Fig. 3. — The pairwise velocity dispersion determined from the redshift power spectrum for 
the Las Campanas Redshift Survey. The whole sample, the south subsample, and the north 
subsample are represented by the triangles (with error bars) , the dashed and the dotted lines 
(without error bars) respectively. 
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Fig. 4. — The real space power spectrum determined from the redshift power spectrum for 
the Las Campanas Redshift Survey. The whole sample, the south subsample, and the north 
subsample are represented by the triangles (with error bars) , the dashed and the dotted lines 
(without error bars) respectively. 
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Fig. 5. — The inverse of the damping function D{k,ji) of the cluster-weight biased 
mock samples with Q = 0.2, plotted as a function of kficr v . In this plot, we have set 
a v = 600 km s -1 . 
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Fig. 6. — The pairwise velocity dispersion and real-space power spectrum determined for 
the 60 mock samples. The error bar is the expected la error of the mean value of the 
60 samples. The solid lines are the pairwise velocity dispersion measured from the peculiar 
velocity (assuming k — 1/r for the plot) and the real-space power spectrum from the particle 
spatial coordinates. Similar to the solid lines, the dashed ones are measured from the full 
simulations but without the cluster-weighted bias. 
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Fig. 7. — The pairwise velocity dispersion and real-space power spectrum determined for the 
Las Campanas Redshift Survey, compared to three CDM models with the cluster-weighted 
bias. The triangles are the results of the LCRS, and the lines are the la upper and lower 
limits of the CDM models. 



